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JTl ■ ABSTRACT 
O 

p We study the effects of dust grain size on the spectral energy distribution 

Q ■ (SED) of spherical circumstellar envelopes. Based on the self-similarity relations 

^ . of dusty SEDs derived by Ivezic & Elitzur (1997), we expect an approximate 

■ invariance of the IR SED for models with different grain sizes. Approximate 

invariance follows from the fact that differently sized grains have similar optical 
^ . properties at long wavelengths where the dust reprocesses the starlight. In this 

. paper, we discuss what are the physical requirements on the model parameters 

to maintain the approximate invariance of the IR SED. Single grain size models 
are studied for a wide range of grain sizes in three optical depth regimes: opti- 
cally thin models, moderate opacities, and very optically thick models. In this 
study, we find limits for the cases where the IR SED is and is not capable of 
conveying information about grain sizes, and to what extent it does so. We find 
that approximate invariance occurs for a much larger range of grain sizes than 
previously believed, and, when approximate invariance holds, the SED is con- 
trolled mainly by one parameter, the reprocessing optical depth, a quantity that 
measures the fraction of starlight that is absorbed by the dust grains. Models 
with a grain size distribution are studied as well. For these models, we find that, 
in many instances, the concept of approximate invariance may be extended from 
the IR SED to all wavelengths. This means that, for a wide range of optical 
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depths, models with different grain size distributions will produce very similar 
SEDs and, hence, the reprocessing optical depth is the only quantity that can be 
unambiguously obtained from the SED. The observational consequences of this 
result are discussed in detail. Finally, in models with a size distribution, the dif- 
ferent grain sizes each have different equilibrium temperatures. The consequences 
of this effect for the model SED are discussed as well. 

Subject headings: radiative transfer — ISM: dust, extinction — stars: circumstel- 
lar matter — methods: numerical 



1. Introduction 

Dust is associated with many astronomical objects, such as stars and galaxies. Dust 
grains typically absorb radiation at short wavelengths, and since the dust grains are usually 
cooler than the radiation sources, they reemit the absorbed radiation at longer wavelengths. 
Consequently, the presence of dust is usually detected by a flux excess at infrared wave- 
lengths. The usual challenge is to use the spectral energy distribution (SED) of the object, 
to extract information about the nature of the underlying luminous source and the properties 
of the surrounding dust; i.e., its spatial distribution, optical depth and intrinsic properties. 
The intrinsic properties of the dust grains are their chemical composition, condensation 
temperature, shape and size. 

Knowledge of the dust grain sizes present in circumstcllar matter is desirable for many 
astrophysical situations. For example, it is believed that about 60% of the dust injected 
into the ISM originates from the winds of AGB stars (Gehrz 1989), so knowledge of the 
size of the dust grains present in the winds of such stars is fundamental for understanding 
the properties of interstellar dust. Similarly, grain size is important for understanding the 
details of grain formation and mass loss mechanisms for cool stars. If we assume that the 
source spectrum and the spatial distribution of the dust are known for these objects, the 
question becomes to what extent is it possible to obtain the intrinsic properties of the dust 
from knowledge of the SED alone. In this paper, we investigate how grain size affects the 
SED and to what extent the SED can be used to constrain the grain sizes. 

It is well known that models of the SED are not unique. For this reason it is impera- 
tive to identify the fundamental parameters controUing the SED. Recently, Ivezic & Ehtzur 
(1997), hereafter IE97, investigated the scaling properties of the radiation transfer problem 
in spherically symmetric dust-envelopes. They found that the radiation transfer problem in 
dusty media depends only on the following quantities: 
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i) The geometry (shape) of the system (i.e., all distances are proportional to a scale 
factor) ; 

ii) The sublimation temperature of the material; 

iii) The spectral shape of the input radiation, XFx/ F; 

iv) The shape of the dust absorption and scattering opacities, kx/ k^^ and cx/crxo, re- 
spectively, where Aq is the fiducial wavelenght; 

v) The dust scattering phase function (SPF)^; 

vi) The overall optical depth at the fiducial wavelength. 

The above list can be regarded as a set of invariance requirem,ents. Two models that 
have different physical parameters, but meet the above requirements exactly, are equivalent 
and have the same SED. For example, the physical dimensions of the dusty region can be 
freely changed (say by changing the stellar luminosity, which increases the dust condensation 
radius) with no effect on the SED, as long as the overall optical depth and shape factors do 
not change. If an invariance requirement is violated, however, then the models are no longer 
equivalent and the SEDs are expected to be different. For instance, if one were to change 
the grain size, the shape of the opacity and the SPF would change, violating requirements 
(iii) and (v). 

For a given class of astronomical objects (e.g., AGB stars), most of the quantities above 
are likely to be similar (geometry and grain composition), so the optical depth becomes the 
single most important parameter that controls the SED. 

It is reasonable to expect that if an invariance requirement is weakly violated, the SED 
will still be approximately the same. This was noted by IE97 who pointed out that if the 
grains are very small (about a tenth of wavelength of the peak of the source spectrum) then 
the shape of the opacities are very similar and the grain size is irrelevant for the problem. 
For a 3000 K source, for example, the upper limit for the grain size is about 0.05/im. Grains 
larger than this upper limit will, according to IE97, significantly alter the results. 

In the case of optically thin envelopes, this condition on the maximum grain size can 
be further relaxed, because the star completely dominates the optical SED while the grains 
produce the IR SED. For the IR SED to remain similar between models with different grain 
radii, it is evident that both the shape of the IR emissivity and the reprocessed luminosity 



^This condition was not required by Ivezic & Elitzur (1997) since they only considered isotropic scattering. 
We consider dust with anisotropic phase function, therefore this additional requirement is necessary. 
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must be the similar. For the shape of the IR emissivity to be the similar, it is necessary 
that the grains have similar absorption efficiencies only in the IR; this relaxed requirement 
increases the maximum grain size to at least 15fim, assuming the grains have equilibrum 
temperatures of about 1000 K. The condition that the reprocessed luminosity be similar is 
equivalent to imposing a condition on the optical depth (see next section). 

This example shows that, under certain circumstances, it is possible to relax the invari- 
ance requirements and still maintain the invariance of the SED, at least approximately. We 
call this concept approximate invariance, and in this paper we demonstrate how it reveals 
the similarities between models with different grain sizes. In particular, wc study to what 
extent the grain sizes can be changed, while maintaining the approximate invariance of the 
SED under a wide range of model parameters. For each case presented, we identify the 
primary parameter (like optical depth) that controls the SED. 

In the next section, we explore in more detail the physical requirements for approximate 
invariance when the dust grain size is varied. In section 3, we briefly describe the Monte 
Carlo code used for the calculations. In section 4, we present the effects of grain size on the 
SED for single grain size models. Section 5 extends the study to models with a grain size 
distribution, and a discussion and summary of the results is presented in section 6. 

2. Approximate Invariance 

Consider an astrophysical system, which consists of a luminous source (e.g., a star) 
surrounded by a dusty envelope. The physical description of this system must include both 
the shape and spectrum of the source and the dust properties (chemical composition, grain 
size, spatial distribution and optical depth). One characteristic of such a model is the 

presence of an inner cavity. Typically, for a star that is losing mass (e.g., an AGB star), 
the location of the inner cavity will be controlled by dust sublimation/condensation. In 
contrast, a star that has stopped losing mass (e.g., a planetary nebulae) may have a much 
larger cavity, with a correspondingly cooler radiative equilibrium temperature at the inner 
edge of the cavity. 

Now consider different models of the system, where we vary only three parameters: the 
dust grain size (a), the cavity inner radius (ri), and the optical depth (r). Such models 
violate some of the invariance requirements of section 1, thus they should have different 
SEDs. However, it follows from the idea of approximate invariance (AI) that their IR SED 
will be similar, provided that the three conditions bellow are satisfied: 

1 . The shape of the absorption efficiency factor for the different grain sizes is similar in 
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the spectral range where the grains emit most of their thermal flux; 

2. The integrated (bolometric) IR luminosity is the same; and 

3. The temperature of the grains at the inner edge of the envelope is the same for all 
grain sizes. 

In the following, we investigate the requirements these conditions impose on the three 
model parameters (optical depth, grain size and cavity radius). Let us first establish the 
spectral region where we expect the absorption efficiencies Qabs of differently sized grains 
to be approximately equal (condition 1). Figure 1 shows Qabs for cosmic silicate of various 
grain sizes (optical data from Ossenkopt, Henning & Mathis (1992)). Mie theory defines 
three spectral regions where the grains have common properties: 

i) A < vra. In this region, called the geometrical limit, the optical properties are roughly 
independent of A, and Q^hs ~ const.; 

ii) A > 5/Ta. In this interval, the efficiencies arc in the so-called Rayleigh limit, where 
the size of the particle is much smaller than the wavelength. In this limit, Qahs oc 1/A for 
amorphous grains and Qabs oc 1/A^ for crystalline grains; 

iii) Tra < A < Stto. In this region, the efficiencies are calculated using Mie theory 
(spherical particles) and are dependent on grain size. We call this interval the Mie region. 

These three regions are easily recognizable in Figure 1. We see that all grains have 
approximately the same efficiencies for short wavelengths (region i) and that the shape of 
the efficiency factors are the same for A > 10 /xm (region ii). In the visible and NIR region, 
the efficiencies are very different. From the above discussion, it is easy to see that condition 1 
of AI is satisfied when most of the IR radiation is emitted in the Rayleigh limit. In practice, 
this sets a constraint on the maximum grain size that satisfies condition 1 of AI. For a given 
temperature of the inner cavity, Tj, the largest grain size satisfying condition 1 of AI is given 

by 

a<lS5/T,pim. (1) 

For example, for Ti = 1000 K (approximately the sublimation temperature for most grain 
species), the maximum grain size satisfying condition 1 of AI is about 0.2 //m. 

Condition 2 of AI requires that the emergent envelope luminosity, Lj-cp, remain fixed. 
The envelope is powered by absorbing some fraction of the stellar luminosity, so condition 
2 is equivalent to requiring that this absorbed fraction remain fixed. To characterize this 
condition, we introduce the quantity Trep, the reprocessing optical depth. This quantity is 
defined so that envelopes with the same Trep will reprocess (i.e., absorb) the same fraction of 
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the input energy. The reprocessing optical depth is thus given by the condition 

1 _ e-^-p = Lrep/L^ . (2) 



To define the quantity Lfcp it is useful conceptually to divide the emergent radiation into 
stellar photons and envelope photons. In a Monte Carlo simulation, the stellar luminosity 
is divided into N equal energy photon packets (note that packets with different frequencies 
contain a different number of physical photons) . The energy per packet is given by 



(3) 



where L^, is the stellar luminosity and At is an arbitrary simulation time. Stellar photons 
are emitted from the star and they propagate through the envelope, where they may be 
scattered or absorbed. If a stellar photon is absorbed, radiative equilibrium requires that 
it be rcemitted as an envelope photon. Note that scattering does not change the photon 
type (stellar vs. envelope). The fraction of the stellar luminosity reprocessed (absorbed) 
by the envelope may easily be determined by counting the number of stellar photon packets 
absorbed by the envelope, Nf*^. This is equivalent to counting the number of stellar photons 
that emerge from the envelope without absorption, N^'^. Thus the luminosity reprocessed 
by the envelope is 



Nf'E^/At ={N- N^"')E^/At = 1^,-11 



^rep 



(4) 



where L^°^ is the emergent stellar luminosity. This is the method we use to measure Lrep 
in this paper. Note that scattering affects L^ep because a photon packet that is scattered 
can be absorbed subsequently. Owing to these scattering effects, it is not straightforward to 
obtain Lj-cp from standard radiative transfer quantities. The reprocessed luminosity is given 
by the emergent luminosity of the envelope, so 



^rep 



p 



(5) 



where da/dD, is the differential scattering cross section, which is a function of both incoming 
and outgoing directions (h' and n, respectively), and /|"^ is the envelope contribution (Monte 
Carlo envelope photons) to the specific intensity (this includes scattered envelope photons 
but docs not include scattered stellar photons). Energy conservation requires that L^^p equals 
the energy absorbed in the envelope, so L^ep can also be written as 



Lrep = 47r y KxpJ^dXdV , 



(6) 



where J\ is the stellar contribution to the mean intensity, including scattered starlight. 
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For low optical depths in the point source approximation (r^ ^ B^), J\ — > where 
H\ is the source flux. Defining H* = H\d\ and the flux mean opacity 

«F - , [J) 

eq. (6) can be rewritten, in the optically thin limit, as 

Lrep ^ 4^ y i^xpH^dX A'Kv'^dr = IGtt^ j K^pH^r'^dr ^ ^* J i^Fpdr = L^Tp , (8) 

where tf is the flux mean optical depth. From eqs. (2) and (8) we see that in the optically 
thin limit, r^cp is equivalent to Tp. 

For larger opacities, I^"^ of eq. 5 is not zero, and its determination must take into 
account the full radiative transfer problem, including multiple scattering effects. In general, 
multiple scattering will increase the effective opacity of the envelope, making r^ep larger than 
Tp. From this it follows that the detailed shape of the scattering phase function (SPF) can 
be of importance for determining the amount of reprocessing in the wind. That is the main 
reason why the SPF is listed above (item v in section 1) as one of the quantities necessary 
to completely specify the radiative transfer problem in dusty media. In section 4.4 we will 
expand on this issue with further details. 

An important consequence of eq. (5) it that, because Lrep depends on the details of the 
radiative transfer, it is not straightforward to obtain a relation between Lrep (or '^rep) ^-nd 
the optical depth at the fiducial wavelength (say, Tv)- In this paper we demonstrate that, 
because of AI, r^cp is the most suitable parameter for SED classification; it follows, then, 
that the optical depth Ty is not a suitable parameter, contrary to the usual approach in the 
literature. 

Although Trep (or Lrep) is a well defined quantity, relating this quantity to the physical 
parameters of the model (such as mass loss rate) is a rather involved problem, requiring the 
full solution of the radiative transfer. On the other hand, in many instances Lrep is a well 
defined observational quantity. This is the case, for example, for optically thin envelopes 
surrounding a hot star, in which the attenuated stellar spectrum can be easily separated from 
the IR emission of the grains. In the following, where we discuss the detailed consequences of 
AI, we will show that this issue is at the heart of the fundamental uncertainties in determining 
mass loss rates using the SED. 

Finally, condition 3 of AI determines the radius of the cavity. In particular, the radius 
will be the dust destruction radius when the presence of the cavity is controlled by dust 
sublimation/condensation. Note that by specifying the cavity radius in stellar radii we are 
not introducing another radial scale in the problem. 
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3. The Monte Carlo Code 

For the calculations shown in this paper, we have used a Monte Carlo code developed to 
solve the general problem of the radiative transfer plus radiative equilibrium in dusty media 
(Carciofi 2001). The code uses a standard Monte Carlo simulation, which basically follows 
the path of a large number of photon packets as they are scattered, absorbed and reemited 
within a prescribed medium (for further details see, e.g.. Code & Whitney (1995)). The 
radiative equilibrium is solved using the method described in Bjorkman & Wood (2001), 
in which the grain temperature and emitted spectrum are corrected in the course of the 
simulation. 

Although the radiative equilibrium calculation does not require any iteration, quantities 
like the dust condensation radius ri and the reprocessing optical depth Trep (see section 2) do 
require iteration, since they depend on the results of the radiation transfer. Each iteration 
returns a temperature of the inner cavity and the amount of reprocessed energy. The values 
of Ti and the optical depth are modified accordingly, and a new simulation is run, until the 
correct values are obtained. 

The code is fully 3-D, so it is capable of solving the radiative transfer problem for 
many circumstellar geometries and density distributions, including the presence of multiple 
sources. The transfer of polarized radiation is included, so the code is capable of providing 
the entire SED, the polarization as a function of wavelength, images at specific wavelengths 
as well as polarization maps. 

The code treats an arbitrary mixture of grain species (amorphous carbon, silicate, etc.) 
and a different size distribution can be assigned to each grain species. The grain scatter- 
ing and absorption properties are calculated using Mie theory, without approximations, to 
obtain the correct form of the scattering phase function for polarized incident radiation. 
The radiative equilibrium condition is imposed separately for each grain type and size, so 
the code calculates independent radiative equilibrium temperatures for each grain type and 
size. This procedure contrasts with the usual procedure in the literature, which consists of 
assigning a single temperature to all grain sizes, whose opacity and scattering properties are 
obtained by averaging over the size distribution. In section 5, we will show that, for hot 
stars mainly, individual grain sizes can have very different equilibrium temperatures, which 
has important effects for the model IR SED. 

In summary, the only approximation used in our dust model is that the shape of dust 
grain is spherical. 
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4. Single Grain Size Models 

In this section, we explore the effects of grain size on the SED for single grain size 
models, and interpret the results in the framework of approximate invariance. Our goal is 
to determine to what extent the grain size can be varied while maintaining the approximate 
invariance of the SED. We begin by defining the basic model we use in this and in the next 
section. It consists of a central star of unit radius that emits a black body spectrum of 
temperature Teflf, surrounded by a spherical dust shell with internal radius ri and external 
radius 1000 ri. The internal radius is set by the temperature of the cavity wall, Ti, often 
taken to be the dust condensation temperature. The dust density profile is proportional 
to r~^, and the dust grains are spherical. This is a simple model, but useful with respect 
to many astrophysical situations, such as the nearly spherical circumstellar envelopes found 
around many RGB and AGB stars. 

Three different chemical compositions were studied: cosmic silicate, with optical con- 
stants given by Ossenkopt, Henning & Mathis (1992), amorphous carbon (optical constants 
by Zubko et al. (1996)) and silicon carbide (optical constants by Pegourie (1988)). The sub- 
limation temperature for these materials was taken to be 1000 K for silicate, 800 K for amor- 
phous carbon and 1500 K for silicon carbide. These numbers are somewhat arbitrary, but 
they reflect approximately the temperatures used in the literature (e.g., IE97, and Lorentz- 
Martins & Lefevre (1994)). 

The model parameters we vary are the dust properties, optical depth, source tempera- 
ture, and temperature of the cavity wall. To comply with the three conditions of AI, stated 
in section 2, we compare models with the same chemical composition, the same temperature 
of the cavity wall, and the same reprocessing optical depth. 

We shall distinguish between three different cases, corresponding to three optical depth 
regimes, which will be studied separately. In case A, the envelope is optically thin, both in 
the wavelengths emitted by the source (the extinction region) and in the IR (the reprocessing 
region). In case B, the envelope is optically thick in the extinction region, but is optically 
thin to reprocessed radiation, and in case C the envelope is optically thick in both regions 
(up to the wavelength for which the Rayleigh limit is achieved for all grain sizes considered) . 

4.1. Case A: Optically Thin Envelopes 

In this case, we determine the V-band optical depth Ty for each model using the con- 
dition Trep = 0.1, which corresponds to reprocessing 9.5% of the input radiation. In the 
extinction region, about 90% of the flux is stellar in origin, so we expect SEDs from different 



-10- 



models will be similar in the extinction region, at least at the 10% level. The IR SED will 
also be similar, as a result of the AI. Hence, we expect to observe a similarity between all 
models for all wavelengths. 

We illustrate this situation in Figure 2. Each panel shows the SED for four different 
models, each with a different grain size, for cosmic silicate grains with a given stellar tem- 
perature. In Figure 3, wc show the extinction optical depth as a function of wavelength for 
the models with T^s = 2500 K. The other model parameters, cavity radius, grain mass and 
V-band optical depth are listed in Table 1. 

Figure 2 corroborates our above qualitative discussion. For all grain sizes and stellar 
temperatures, the short wavelength side of the curve is very similar, as a consequence of the 
low optical depth combined with the same total attenuation. Similarly, the overall amplitude 
of the IR emission scales approximately with Tj-cp and the shape is set by the grain absorption 
efficiency, which is the same for each grain size. As a result, the overall SED is very similar 
for all models. The mean difference between the models with Tes — 2500 K is about 7%, and 
for the models with Tgs — 20000 K, the mean difference is 9%. The similarity between the 
model SEDs, not only in the IR but also in the extinction region, allows us to extend the 
consequence of AI. We conclude that in the optically thin limit the entire SED is similar 
when the conditions of AI are met. 

The conditions for approximate invariance are violated when we consider grains larger 
than about a — 0.25 iim. These grain sizes fail to satisfy condition 1 of AI, so large differences 
in the IR SED of these models are observed. This is shown in Figure 4, where we plot the 
SEDs for models with grain sizes 0.005, 0.50 and 1.0 //m. The large differences in the IR SED 
are apparent. However, as shown in eq. (1), the largest grain that satisfies condition 2 of AI 
depends on the temperature of the cavity wall. Consequently, if we lower the condensation 
temperature, these larger grains will also satisfy the AI conditions. This is illustrated in 
Figure 5, which shows the SED for envelope models with Ti = 300 K. These plots show that 
a distinction can hardly be made between envelope models with grain sizes now ranging up 
to 1 iim. For T^s — 2500K the mean difference between the SEDs is at most 3%. This has 
the consequence that the SED is a poor constraint for grain sizes, at least for sizes that 
satisfy the conditions of AI. 

We conclude that, for case A, given a grain composition, a temperature for the inner 
cavity, and a spatial distribution of the dust, the SED depends on a single parameter, Trep. 
Another important result is that Trep is not uniquely related to the F-band optical depth 

(or any other wavelength), as shown in Figure 3. All curves are for Trep = 0.1, and yet show 
very large differences in the optical depth for most wavelengths. From this we can conclude 
that the optical depth is not a suitable parameter for SED classification; if we were to group 
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the case A models according to ry, for example, the models would have large differences in 
the IR SED, leading us to the wrong conclusion that the grain size does have an important 
effect on the SED, which we have demonstrated not to be the case. 

The optical depth depends on the product of the grain opacity (which depends on the 
grain size) and the grain density. The fact that grain sizes are very difficult to determine 
from the SED implies a corresponding uncertainty in the dust mass of the envelope, as shown 
in table 1. This uncertainty directly affects the ability to measure the mass loss rates from 
stars using the SED. 

4.2. Ccise B: Optically Thick Extinction with Optically Thin IR Emission 

In this case, the envelope is optically thick in the UV and visible, and optically thin in 
the IR. Based on AI, we expect the behavior of the IR SED to be similar to case A; i.e., the 
shape and level of the IR emission should be the same for all models with the same Tj-^p. In 
contrast to case A, the larger optical depths result in a much larger attenuation of the stellar 
flux. Therefore, the shape of the SED in the extinction region is controlled by the shape of 
the extinction optical depth, which strongly depends on grain size. 

In Figure 6 we present the SED for the same envelope models as in Figure 2, but we 
maintain Xrcp = 1. This implies 63.2% of the input radiation is reprocessed by the envelope. 
As before, other models parameters are listed in Table 2, and in Figure 7 we display the 
extinction optical depth as a function of wavelength. 

We see in Figure 6 that indeed the IR SED for the different models arc similar, while the 
SEDs for short wavelengths arc very different because they depend on the shape of the grain 
efficiencies. It follows that case B is different from case A, in the sense that the SED does 
depend strongly on the dust grain size. There are, however, two points to be considered. 
First, in order to obtain information about grain size, one has to observe the SED in the 
wavelengths where the effects of grain sizes are important (i.e., UV, visible or NIR, depending 
on the model). Second, single grain size models are very crude approximations for actual 
stellar envelopes, where a distribution of grain sizes is expected. In section 5 we study this 
case and show how it affects our conclusions for case B. 

At this point it is useful to compare our results with those of previous authors. For 
example, IE97 show SEDs in their Figure 10 for spherical envelopes of silicate and amorphous 
carbon dust grains, of varying sizes and optical depths. In their plots on the right, where 
the results for silicate grains are displayed, the upper two plots can be compared roughly 
with our Figures 2 and 6. In Figure 10 of IE97, large differences between the models with 



a — 0.05 and 0.1//m are evident. By looking at these plots, the reader is led to believe 
that the IR SED for grains with sizes 0.05 and O.lfim is intrinsically different. However, our 
results show that this is not the case. The IR SED of single grain size models whose sizes 
comply with the conditions of AI have very similar shapes and levels. 

The differences between Figure 10 of IE97 and our Figures 2 and 6 owe to the choice 
of the model optical depth. In IE97, all models have the same extinction K-band optical 
depth. Our models, on the contrary, have all the same r^ep but very different ry, as seen in 
Figures 3 and 7. Thus in case B, we also find that Tj-ep is the best parameter to reveal the 
essential similarities of models with different grain sizes (which must be lower than the limit 
set by condition 2 of AI). 



In this case, the envelope is so optically thick that 100% of the input radiation is 
reprocessed by the dust grains. As a consequence, the quantity Tj-cp becomes ill-defined, so 
another parameter must be found to reveal the similarities resulting from AI. 

By definition, all models, in this case, are optically thick to the reprocessed radiation. 
This means that most of the radiation will emerge at wavelengths larger than the Wien peak 
of the cavity wall (i.e., A > 3 /im for Ti = lOOOK). For a given wavelength the envelope can 
be divided into two different regions with different properties concerning the radiation field: 
an inner region with radius between ri and rr=i, and an outer region, between rr=i and r^, 
where re is the envelope outer radius, and rT-=i is given by 



i.e., r^=i is the point where the radial optical depth of the envelope is unity (this quantity 
depends on the wavelength). By definition, radiation emitted inside the inner region is 
unlikely to leave the envelope because the optical depth is larger than 1; hence, a reasonable 
approximation for the emerging IR fiux is given by the volume integral of the radiation 
emitted by the outer (optically thin) region 



In order to find an analytic expression for Fx, we will suppose that both the temperature 
in the outer part of the envelope and the grain opacity can be approximated by power-laws, 
so we write 



4.3. Case C: Optically Thick Envelopes also in the IR 




(9) 




(10) 




(11) 
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and 

Kx^Ko(jJ , (12) 

where Aq is an arbitrary constant. The assumption for the opacity is very reasonable, because 
by definition most of the flux emerges at wavelengths larger than the Rayleigh limit (see 
Figure 1 and discussion in section 2). The assumption for the temperature is justified by 
the fact that the outer part of the envelope is optically thin; it can be shown that, if the 
dust opacity is given by eq. (12), the radial dependence of the temperature of an optically 
thin envelope is T(r) oc r"^/^"*"^) (e.g. Martin (1978), eq. 7.3). Our detailed calculations, 
shown below, validate this assumption. Finally, assuming a dust density p(r) , also given by 
a power-law 

/r - \ n 

P(r)=Po(^) , (13) 
we calculate r^=i using eqs. (9), (12) and (13), which gives 



n V A 

where 



n-l 



ro-\ (14) 



is the envelope optical depth at A = Aq. 

Substituting eqs. (11), (12), (13), and (14) in eq. (10) and letting Vg — > oo, we obtain 
the spectral shape of the radiation emitted by the outer part of the envelope 

^-KX"^---^ T-^du, (16) 



where 

n— 1 + sp 

u, = rl-) , (17, 



and 

The normalization constant K is defined by the condition 



r^T^vf^- (18) 



\F,d\ _ ^ 



In eq. (16) the only parameter that depends on the grain size is F in the lower integration 
limit, Ui. It follows that, if the values of Tq and Tq of models with different grain sizes are 
such that F is the same, these models will produce the same SED. 
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We show in Figures 8 and 9 the SEDs and grain equihbrium temperatures for 20000 K 
models with different grain sizes. Figure 9 shows that the temperature power-law exponent 
has a value s = 0.4 at large r. Other model parameters arc listed in Table 3. Note that the 
optical depths are very high (ranging from Tq = 3.4 to 5.0 for Aq = 100 /im), according to the 
assumption for case C. These optical depths were chosen so that the parameter F is the same 
among the different models (i.e., we measured Tq and adjusted tq so that TqTo"" = constant). 
The resulting SEDs are very similar for A > 20/im, which indicates that the parameter F 
is an adequate scaling parameter for very high optical depths. For shorter wavelengths the 
SED invariance breaks down, because the opacities are no longer described by a power-law, 
and the SEDs are different. 

Also shown in Figure 8 is the analytic expression for the emerging flux, eq. (16). This 
expression reproduces well the shape of the SED for long wavelengths; for shorter wavelengths 
the agreement is not good for the reason stated above: our assumption for the grain opacity 
fails in the well-known silicate spectral features at 10 and 20 /im, and as we go into the 
optical limit. 

We conclude that, given a dust grain composition and spatial distribution, the shape 
of the IR SED is controlled by a single parameter, F. From the observational point of view, 
wc have a situation similar to case B, where the IR emission bump does not convey any 
information about the dust grain size. However at shorter wavelengths, AI breaks down, so 
one can determine information about the dust grain size, if there is sufficient observable flux 
at these wavelengths. 



4.4. Effects of the Scattering Phcise Function on the SED 

In the previous sections we studied in detail the effects of grain size on the SED, and 
determined to what extent the grain size (i.e., opacity) can be changed while maintaining 
the approximate invariance of the SED. An additional effect is that, when the grain sizes 
are changed, both the opacity and the SPF change. Changing the SPF breaks invariance 
requirement (v) of section 1 and thus will alter the SED. The results shown in the previous 
sections were calculated using Mie theory which has an anisotropic SPF, but a common 
simplification is to use an isotropic SPF. In this section, we briefly discuss how an isotropic 
SPF affects our previous results. 

In general, the primary change is the l/-band optical depth required to reproduce Trep; 
the SED is largely unaffected (except at the few percent level). In Table 4, we compare the 
y-band optical depth between models with an anisotropic SPF (r^) and an isotropic SPF 
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{Ty°) for both case A and B. One should compare column 1 with column 2, and column 4 
with column 5. Because all models are subject to the condition Trep = 0.1 (case A) or 1 
(case B), the difference between Ty and ry° for each case gauges the effects of the SPF on 
the amount of reprocessing for each model. 

For case A models, Ty and Ty° are roughly the same (at the 10% level). This is to be 
expected because, in the optically thin limit, 7^°"^ in eq. (5) is small, so the SPF (da/dQ) 
has little effect on Lrep. 

For case B models, the differences between Ty and Ty° are larger, but still not very large 
(at most 22%). We conclude that for optically thin models the SPF is largely irrelevant, 
while for optically thick models, the SPF affects the conversion from Ty to Trep, but at most 
to a level at 20% to 30%. 

Also shown in Table 4 is the flux mean optical depth, tf, deflned in eq. 8. One should 
compare these values with the corresponding value of Trep. The difference between tf and 
Trep arises from multiple scattering effects, and for this reason they are much smaller for 
optically thin models than for optically thick models. Many authors use Tp as a parameter 
for SED classification; the results in Table 4 show that this is not a good parameter, at least 
for optical depths in the range of case B models, for which multiple scattering effects are 
important. 



5. Models with a Grain Size Distribution 

The single grain size assumption can be a poor approximation for describing the dust 
of an astrophysical system. For example, fitting the interstellar extinction curve requires 
a distribution of grain sizes. Mathis, Rumpl & Nordsieck (1977, hereafter MRN) used a 
simple power law size distribution with index q — —3.5, and a lower and upper cutoff size, 
Omin ~ 0.05 /im and amax ~ 0.25 /im. This distribution has been revised many times in the 
hterature. For example, Kim, Martin & Hendry (1994) proposed a distribution in the form 

fin 

— oca-«e-«/"%a>ami„ , (20) 

where Oq = 0.14 //m for silicate grains and Oq = 0.28 //m for graphite grains. The upper 
cutoff was replaced by an exponential decay, to allow for the presence of large dust grains 
{a > 0.2 urn). 

Although the fact that the dust has a distribution of grain sizes is well established, the 
exact form of this distribution is still an open issue. As we will demonstrate below, part 
of this uncertainty is due to the AI of the problem, which makes two models with similar 
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size distributions virtually indistinguishable. Another source of uncertainty lies on our yet 
incomplete knowledge of the details of grain formation and growth. To complicate matters 
even further, the grain size distribution is likely to vary across the wind because different 
grain sizes have different drift velocities (Elitzur & Ivezic 2001). 

In this section, we study the SED of models with a grain size distribution. For simplic- 
ity, we adopt a power-law (MRN) grain size distribution function; although the quantitative 
details of the results shown below do depend on the adopted size distribution, the funda- 
mental concepts discussed do not. Previously, in section 4, we used the concept of Al to 
reveal important similarities in the IR SED (and in the extinction SED for optically thin 
models) for single grain size models. In the following, we study what controls the spectral 
shape of the SED for models with a size distribution in an attempt to identify the same sort 
of similarities. 



Let us consider an envelope model, consisting of spherical dust grains of the same 
chemical composition and sizes distributed according to a given size distribution function, 
dn/da. The envelope IR SED, Fir, will be given by the sum of the contributions from each 
grain size 



where /(a) is the wavelength integrated emission for each grain size, a measure of how the 
absorbed energy is split between the different grain sizes, and J-\-^{\, a) is the spectral shape 
of the spectrum emitted by the grains with radius a. 

The simplest situation is when the spectrum emitted by all grains has approximately 
the same spectral shape (i.e., J?-ir(A, a) is independent of a). In this case, the spectral shape 
of the IR SED of the envelope will be similar to the spectral shape of the individual grains, 
regardless of how the absorbed energy is split between the different grains. A somewhat 
similar situation occurs when a particular range of grain sizes, with a similar emission spec- 
trum, completely dominates the emission. Here, even if the spectral shape of the other grain 
sizes is different, it does not affect the envelope IR SED. An important difference from the 
first situation is that, here, the manner in which the absorbed energy is distributed between 
the grain sizes is important. For both these cases, we expect a close similarity to the corre- 
sponding single-sized grain model. An opposite situation occurs when different grain sizes 
emit a different spectral shape with similar integrated emissions. In this case, the envelope 
IR SED can no longer be described using a single grain size model. 



5.1. Spectral Shape of the SED 




(21) 
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What controls ^ir(A, a) is the spectral shape of the absorption efficiencies and the grain 
temperature. From the conditions of AI, it is evident that the grain sizes with similar emission 
spectra will be the ones that comply with conditions 1 and 3. The first condition requires 
all grain sizes of the distribution must be smaller than a maximum size, given by eq. (1). 
The third condition implies that all grain sizes must have similar equilibrium temperatures. 
However, owing to their different absorption efficiencies, different grain sizes can have very 
different equilibrium temperatures. This effect can be inferred from the second column of 
Table 1: different grain sizes have different condensation radii, which indicates that, if these 
grains were to coexist at the same point in space, the grain sizes with larger condensation 
radii will be hotter than those with lower condensation radii. Recently, Wolf (2003) studied 
the condensation temperature of the individual grain sizes in a mixture of different grains 
sizes, both in 1-D and 2-D dust shells, and found that the temperature difference spans a 
range of up to ~ 250 K, although this value is highly dependent of the choice of the dust 
properties. 

It is easy to show that, in the Monte Carlo radiative equilibrium scheme, the temperature 
T{a) of the grains with radius a at a given point of the envelope will depend on the number 
of photons absorbed by these grains, iVa,bs(a), and on the Planck mean opacity, Kp{T, a) (see 
Bjorkman & Wood (2001), eq. 



T{ay (X 



4 ^abs(a) 



Kp{T, a)' 



The number of absorbed photons is proportional to the absorption cross section averaged 
over the incident photons, so the equation above can be rewritten as 



na)' oc ^ , (22) 
Kp{T, a) 

where 

(.(a)) ^ LJ^MM^ . (23) 

The ratio Kx{a)/Kp is shown in Figure 10 for silicate grains with sizes ranging from 
a = 0.005 to 1 /xm, where the Planck mean opacity was calculated for T = 1000 K, the 
condensation temperature of silicate grains. This plot indicates the relative equilibrium 
temperature for different grain sizes as a function of the wavelength of the illuminating 
radiation field. For example, if the grains are illuminated by a black body radiation of 
temperature T^s — 2500 K, which peaks at A = 1.2 /xm, the grains with intermediate sizes 
(a 0.25 to 0.50 //m) will be the hottest. Note, however, the range of temperatures will be 
relatively small in this case. On the other hand, if the grains are illuminated by a 20000K 
black body radiation field, the hottest grains will be the ones with a 0.05 /xm and all the 
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other grain sizes will have much lower equilibrium temperatures, including the very small 
grains (a = 0.005 /xm). Prom eq. (22) and Figure 10, the ratio between the maximum and 
minimum equilibrium temperature can be estimated to be Tmax/^min ~ 2.5. Hence, if the 
hottest grains have a temperature of 1000 K, the coolest ones will be as cold as 400 K. 

This temperature difference between the different grain sizes is to be expected only 
for optically thin models. For optically thicker models, the temperature differences will be 
smaller, and for very optically thick envelopes the temperatures will be the same, regardless 
of the grain size. This occurs because at very large optical depths, Jx = Bx, and the grains 
arc in thermal equilibrium with the radiation field. 

We now can see an important difference between models of cool and hot stars. For cool 
stars, the equilibrium temperatures of the different grain sizes are not very different, so we 
expect the shape of the IR SED to be similar to those of the individual grain sizes. For 
hot stars, however, the spectral shape for each grain size will be very different. From the 
previous discussion, we infer that the IR SEDs of hot stars will depend on how the flux is 
split among the different grain sizes. 

For an optically thin model, it is possible to estimate the values of /(a), the wavelength 
integrated emission for grains of size a. From Kirchhoff 's law, the emission must be the same 
as the wavelength integrated absorption, i.e.. 



This relative emission is plotted in Figure 11 for a standard MRN size distribution func- 
tion, with q = —3.5. This figure shows the relative /(a) as a function of the illuminating 
wavelength for the same grain sizes shown in Figure 10. Note that, with the exception of 
the smallest grain considered, the curves of Figure 11 scale approximately with the curves 
of Figure 10, which means that, in general, the coolest grains will have the lowest integrated 
emission and vice-versa. As we will see later, this has an important consequence for the hot 
star case, for which the equilibrium temperature of the grains can be very different; however, 
because the hotter grains dominate the emission spectrum, the envelope IR SED will still be 
similar to the spectral shape of the individual grains. 

The smallest grains are an exception for this rule. If we compare them with the grains 
with a ~ 0.05 /im, we see that they reprocess about the same energy, but have lower equilib- 
rium temperatures. This creates a situation in which grains with two different equilibrium 
temperatures both contribute equally to the final emission spectrum. 

The above discussion can be summarized as follows. For cool stars, we expect the IR 
SED of the distribution to be similar to that of the single grain models because all the grains 




(24) 
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have similar temperatures and, hence, similar emission spectra. For hot stars, if we exclude 
the smallest grains (~ 0.005 /xm), the IR SED of the distribution will be similar to the SED 
of the hottest grain, which dominates the emission. It is important to note, however, that 
these conclusions are valid only for the size distribution function assumed above. A different 
value of q, for example, will lead to changes in the relative values for /(a). 



5.2. Approximate Invciriance for Grain Size Distributions 

In this section, wc follow the same approach as section 4 and compare model SEDs 
for different grain size distributions. More specifically, we employ a MRN size distribution 
function with q = —3.5, flmin = 0.05 fim, and different values of Omax- As before, we compare 
models with the same r^cp, to ensure that condition 2 of AI is satisfied. As was seen in 
section 4, different grain sizes can have different condensation radii. However, large grains 
typically form from smaller grains, so, as an approximation, we set the cavity radius to be the 
condensation radius of the smallest grain. In section 5.3, we address the consequences and 
validity of this approach. Finally, we defer discussion of the smallest grains (a ^ 0.005 iim) 
until later, for the reasons mentioned in the last section. 



5.2.1. Case A: Optically Thin Envelopes 

The case A results for models with a grain size distribution are essentially the same 
as case A for single grain size models. If we compare model SEDs for models with a^ax ^ 
0.25 /im, the maximum grain size to satisfy condition 1 of AI, the curves are identical at the 
3 to 5% level, depending on the stellar temperature, as shown in Figure 12. This reinforces 
the conclusion for case A in section 4.1, where we saw that the SED is a poor indicator of 
grain size (as long as the grain sizes present in the distribution each satisfy the condition of 
AI). 

It is interesting to note that, this time, the smallest differences are observed for the 
hot star models, in contrast with the single grain size case, where the smallest differences 
were observed for cool star models. The reason for this can be understood from Figure 11, 
which shows that, for hot stars, a single grain size (a ^ 0.05 /xm) dominates the emission. In 
contrast, all grain sizes contribute to the cool star IR SED with slightly different emission 
spectra owing to small differences in the equilibrium temperatures. Another interesting point 
is that, if we include distributions with Omax larger than that satisfying condition 1 of AI, we 
still obtain very similar SEDs. For example, when Cmax = 0.50 /im, the differences between 
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the model SEDs is less than 10%. This, again, contrasts with the results for single grain size 
models. 



5.2.2. Case B: Optically Thick Extinction with Optically Thin IR Emission 

Figure 13 shows the results for models with Omax = 0.15 and 0.25 //m and optical depth 
set so that Trep = 1. Note the striking similarity in the extinction region, in contrast with the 
previous results for single grain size models (see Figure 6) . The mean difference between the 

two SEDs is now only 11% for Tcff = 2500 K and 10% for Teg = 20000 K. The much closer 
similarity in the extinction region can be understood in terms of the extinction optical depth, 
shown in Figure 14. As a result of the averaging of the optical properties with respect to the 
distribution function, the shape of the extinction optical depth of the two size distribution 
models is very similar. 

The overall similarity of the SEDs for case B has important theoretical and observational 
consequences. From the theoretical point of view, this allows us to extend the consequences 
of AI from the IR SED to the entire spectrum, as was done for case A, but this time for a 
much broader range of optical depths. From the observational point of view, this extends 
the conclusion that SED is a poor indicator of grain size, to include all optical depths. 

5.2.3. Case C: Optically Thick Envelopes also in the IR 

This case, in all aspects, is similar to case C for single grain size models. Since the optical 
depth is very high, the grains tend to thermalize, approaching local thermal equilibrium, 
where no temperature difference is expected for the different grain sizes. Hence, the set of 
different grain sizes will behave as a single-sized grain whose optical properties are given by 
the average of the properties with respect to the distribution function. 

5.3. Brecikdown of Approximate Invar iance 

The similarity of the SEDs for models with different size distributions for all optical 
depth regimes, illustrated in the last section, changes when we consider very small grains. 

Figure 15 shows the case B SEDs for models with amin = 0.005 /xm and Omax = 0.05, 0.15 
and 0.25 /im. As before, the models with large flmax have very similar SEDs in all spectral 
regions (compare Figure 15 with Figure 13); however, a significant difference in the extinction 
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region arises for the model with amax = 0.05 //m. The reason for this difference is easily 
understood from the different shape of the extinction optical depth, shown in Figure 16. 
This difference results from the fact that the smallest grains (in contrast with the larger 
ones) did not reach the geometrical limit throughout the extinction region. 

5.4. Breakdown of Grain Condensation Radius 

In Figure 17 we compare the SED for models with the same a^aax (0.25 //m) and different 
flmin (0.005 and 0.05 /xm) for a stellar temperature of 20000 K. A large difference in the IR 
SED of the two models is observed. This difference arises from our choice of the cavity 
radius as the condensation radius of the smaller grains. For the model with a = 0.05 to 
0.25 fim this choice is consistent, because the grains with a ^ 0.05 fim arc the hottest and 
they dominate the IR emission (see section 5.1). For the model with a = 0.005 to 0.25 /xm, 
however, the smallest grain are cooler than the ones with a 0.05 //m. From the choice of 
the inner radius, it follows that the equilibrium temperatures of the grains with a ~ 0.05 /im 
is about 1250 K, much higher than the grain condensation temperature, which is physically 
impossible. This causes the IR emission to shift to lower wavelengths, as seen in Figure 17. 
It is evident that in this case a more consistent model should be considered, in which the 
intermediate (and larger) grains are not allowed to condense until their radiative equilibrium 
temperature drops below the condensation temperature. 

Strictly speaking, this more consistent treatment should be used for all situations stud- 
ied above, but it would demand a much more complicated procedure and would involve 
additional model parameters. The discussion of section 5.1 helps us set useful limits to when 
the approximation of identical condensation radii is valid for all grain sizes. It will be approx- 
imately valid for cool stars, because the differences between the equilibrium temperatures of 
the grains are not large. It will be valid also for hot stars, when the smallest grain considered 
has the largest condensation temperature. Finally, the assumption fails for hot stars, when 
the smallest grain is cooler and has large integrated emission, comparable with hot grains, 
as seen above. 

We can summarize this as follows: if all grain sizes meet the AI requirements, as is 
generally the case for cool stars, then the approximation of same condensation radius for all 
grain sizes is valid. If the grains do not meet the requirements, but a small range of grain 
sizes completely dominates the IR emission, then the approximate invariance of the SED is 
still maintained. The more consistent treatment is required only when the grain sizes do 
not meet the Al requirements and different grain sizes have similar integrated emissivities. 
However, it is likely when one uses the correct condensation radius for each grain size that 
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even in this last case AI will be recovered. 



6. Summciry 

We have studied the effects of grain size on the SED of circumstellar envelopes with 
dust. To do so we introduced the concept of approximate invariance, as a very useful tool 
for revealing the essential similarities of the problem and for systematically exploring a large 
grid of model parameters. 

The concept of AI follows from the fact that the optical properties of differently sized 
grains is similar in certain spectral regions, and from the idea that if a requirement for SED 
invariance is weakly violated (in this case, the shape of the grain opacity), the SED should 
still be approximately the same. We studied separately single grain size models and models 
with a grain size distribution. For both situations, we studied three optical depths regimes: 
optically thin models {case A), optically thick models in the extinction region [case B) and 
very optically thick models {case C), for which most of the radiation emerges in the MIR to 
FIR. 

Our results for case A show that, given a grain composition and dust density distri- 
bution, the SED of models that comply with the conditions of AI are similar not only in 
the IR, but also in the extinction region. The spectral shape of the SED is controlled by a 
single parameter, the reprocessing optical depth, r^ep, defined in eq. (2). This parameter is 
independent of grain size. Our results for case C are similar but, as shown in section 4.3, 
another parameter (F, defined in eq. (18)) controls the shape of the SED. 

The results for case B differ for single grain size models and models with a size distribu- 
tion. For single grain size models, the SEDs are similar only in the IR, while large differences 
are observed in the extinction region, where the grain opacities depend very much on grain 
size. Models with a grain size distribution, on the other hand, display a striking similarity 
in the extinction region, provided the lower limit of the grain size distribution, Omin, is in the 
geometrical limit in the extinction region. It follows that, if Omin is chosen accordingly, the 
case B SEDs for models with a grain size distribution is controlled solely by the parameter 
Trep (the same result as case A). 

The fact that Trep (or T for case C) is the appropriate parameter for SED classification 
indicates that the usual approach in the literature - using ry or tf - will not reveal the 
general invariance of the SED discussed in this paper. We conclude therefore that the SED 
is generally a very poor constraint of grain size. In most situations, observations must 
resolve differences of the order of a few percent to extract information about the grain size. 
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If we add to this the fact that other model parameters, such as dust spatial distribution, 
composition and optical properties are somewhat uncertain, it becomes apparent that the 
task of extracting information about grain sizes from the SED alone may prove very difficult 
or even impossible for spherical geometries. 

If follows from AI that in most circumstances Trep is the only parameter related to 

the grain opacity that can be unambiguously extracted from, the observations. Using Tj-^p 
to determine dust column density or mass is directly subject to the uncertainty in the 
determination of the grain sizes, which directly affects our ability to measure mass loss rates 
using the SED. 

An important physical effect discussed in section 5.1 is that for models with a grain size 
distribution, the different grain sizes have different radiative equilibrium temperatures. This 
implies that a consistent model should include the correct condensation radius for each grain 
size of the distribution. However, we show in section 5.1 that the approximation of using 
identical condensation radii is reasonable when either all grain sizes have similar condensation 
temperatures or when a particular range of grain sizes dominate the IR emission. The only 
situation where the consistent treatment (each grain size having its correct condensation 
radius) is required is when two or more ranges of grain sizes have different equilibrium 
temperatures and similar integrated emission. This situation was found only for models of 
hot stars with very small grains (a < 0.01 fira) present in the distribution. It is important 
to notice that the results presented here depend on the particular choice of the distribution 
function, because this will control the relative contribution of each grain size to the emission. 
We suggest that a study of a given system could start with an analysis of the chosen dust 
size distribution along the lines described in section 5.1, Figures 10 and 11. Such analysis 
can be very useful because one would know beforehand if the approximation of identical 
condensation radii for all grain sizes of the distribution is a valid one. 

We have presented, here, results only for silicate grains. However, this study was carried 
out for amorphous carbon and silicon carbide grains as well. The results for these materials 
are not shown because, although the specific details of the results are different (see Carciofi 
(2001)), the fundamental concepts discussed here remain vahd. 

A.C.C. and A.M.M. acknowledge support from Sao Paulo State Funding Agency (grants 
98/02238-2 and 01/12589-1). A.C.C. and J.E.B. acknowledge support from NSF grant AST- 
9819928. A.M.M. acknowledges partial support from CNPq. 

The authors wish to thank the referee for his/her very useful comments. 



-24- 



REFERENCES 

Bjorkman, J.E. & Wood, K. 2001, ApJ, 554, 615 

Carciofi, A.C. 2001, PhD Thesis, Universidade de Sao Paulo, Sao Paulo, Brazil 
Code, A.D. & Whitney, B.A. 1995, ApJ, 441, 400 
Elitzur, M. & Ivezic, Z. 2001, MNRAS, 327, 403 

Gehrz, R.D. 1989, in lAU Symp. 135: "Interstellar Dust", ed. L.J. AUamandola & A.G.G.M. 
Tielens (Dordrecht: Kluwer), 445. 

Ivezic, Z. & Elitzur, M. 1995, ApJ, 445, 415 

Ivezic, Z. & Elitzur, M. 1997, MNRAS, 287, 799 (IE97) 

Kim, S.H, Martin, P.O. & Hendry, P.D. 1994, ApJ, 422, 164 

Lorentz-Martins, S. and Lefevre, J. 1994, A&A, 291, 831 

Martin, P.G. 1978, Cosmic Dust: Its Impact on Astronomy (Oxford: Clarendon) 
Mathis, J.S., Rumpl, W. & Nordsieck, K.H. 1977,ApJ, 217, 425 
Ossenkopt, V., Henning, Th. & Mathis, J.S. 1992, A&A, 267,567 
Pegourie, B. 1988, A&A, 194, 335 
Wolf, S. 2003, ApJ, 582, 859 

Zubko, V.G., Mennella, V., Colangeh, L. e Bussoletti, E., MNRAS, 282, 1321 



This preprint was prepared with the AAS IM^jX macros v5.0. 




0.01 0.10 1.00 10.00 100.00 1000.00 

X{fMm) 

Fig. 1. — Absorption efficiency factors for spherical cosmic silicate grains. Each line shows 
Qabs for a given grain radius, as indicated. Optical data for the warm version of Ossenkopt, 
Henning & Mathis (1992). 
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Fig. 2. — Single grain size models for case A, with Tj = 1000 K. Shown are the results for 
four envelope models of cosmic silicate grains, with a = 0.005,0.05,0.15 and 0.25 /.tm. The 
left and right panels show the results for a given stellar temperature, as indicated. For each 
model, the optical depth was adjusted so that Trep = 0.1 (see Table 1). Also shown is the 
stellar spectrum {dotted line). 
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Fig. 3. — Extinction optical depth. The optical depths tx vs. wavelength is shown for the 
models of Figure 2, for the case T^s — 2500 K. Note that tx<1 for the wavelengths emitted 
near the peak of the source flux (1 /xm < A < 3 //m), as required for case A. 
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Table 1. Inner cavity radius, grain masses and F-band optical depth for models of 

Figure 2. 



2500 K 



0.005 


4.5 


1.00 


0.308 


0.05 


4.64 


0.97 


0.640 


0.15 


5.62 


0.85 


2.98 


0.25 


6.26 


0.74 


1.54 


20000 K 


0.005 


1.00 ■ 10^ 


1.00 


2.32 ■ 10-2 


0.05 


1.53- 10^ 


1.00 


2.25 ■ 10-2 


0.15 


1.12 ■ 10^ 


0.92 


0.300 


0.25 


8.12 • 10^ 


0.78 


0.358 



*M/Ma=o.oo5/im is the dust mass normalized to 
the dust mass of the model with a — 0.005 iim 
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Fig. 4. — Single grain size models for case A, with T\ = 1000 K. Shown are the results for 
three envelope models of cosmic silicate grains, with a = 0.005, 0.50 and 1.0 fim. The left 
and right panels show the results for a given stellar temperature, as indicated. For each 
model, the optical depth was adjusted so that r^ep = 0.1. Also shown is the stellar spectrum 
{dotted line). 
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Fig. 5. — Single grain size models for case A, with T[ = 300 K. Shown are the results for six 
envelope models of cosmic silicate grains, with a = 0.005, 0.05, 0.15, 0.25, 0.50 and 1 fira. The 
left and right panels show the results for a given stellar temperature, as indicated. For each 
model, the optical depth was adjusted so that r^ep = 0.1. Also shown is the stellar spectrum 
{dotted line). 
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Fig. 6. — Single grain size models for case B, with Ti = 1000 K. Shown are the results for 
four envelope models of cosmic sihcate grains, with a — 0.005, 0.05, 0.15, 0.25 iim. The left 
and right panels show the results for a given stellar temperature, as indicated. For each 
model, the optical depth was adjusted so that Trep = 1. Also shown is the stellar spectrum 
{dotted line). 
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Fig. 7. — Extinction optical depth. The optical depth tx vs. wavelength is shown for the 
models of Figure 6, for the case with stellar temperature of 2500 K. 
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Table 2. Inner cavity radius, grain masses and V-band optical depth for models of 

Figure 6. 



a(/xm) ri{R^) M/M„=o.oo5 Ty 



2500 K 



0.005 


4.9 


1.00 


3.49 


0.05 


5.3 


1.09 


7.42 


0.15 


7.7 


1.48 


37.0 


0.25 


9.2 


1.30 


17.0 


20000 K 


0.005 


1.04 • 10^ 


1.00 


0.273 


0.05 


1.70 ■ 10^ 


0.96 


0.224 


0.15 


1.17- lO-'^ 


0.74 


2.81 


0.25 


8.61 • 10^ 


0.65 


3.41 
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Fig. 8. — SEDs for five envelope models of cosmic silicate grains, with a = 
0.005,0.05,0.15,0.25 and 0.5 fim. The stellar temperature is 20000 K. Model optical depths 
are given in Table 3. The temperature of the inner cavity was set to 1000 K. Also shown is 
the curve for eq. (16) {thick dotted line). 




Fig. 9. — Temperature profiles for the same models of Fig. 8. Also shown is a fit for the 
temperature of the model with a — 0.05 /im in the outer part of the envelope, with a power- 
law of index s = 0.4. 
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Table 3. Inner cavity radius, grain masses and optical depth at 100 //m for the models 

shown in Figure 8. 



a (//m) 




M/M„=o.005Mm 


TlOO 


0.005 


1.75-103 


1.00 


5.2 


0.05 


1.77-103 


0.99 


5.1 


0.15 


1.85 - 103 


1.03 


4.8 


0.25 


1.97-103 


1.04 


4.3 


0.50 


2.4 - 103 


1.37 


3.8 
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Table 4. F-band optical depth for anisotropic SPF (tv) and isotropic SPF (Ty°), for case 
A (Trep = 0.1) and case B (rrep = 1) models. Also shown is the flux mean optical depth (rp) 

for the models with anisotropic SPF. 



Tv (case A) Ty° (case A) rp (case A) Ty (case B) Ty° (case B) rp (case B) 



2500 K 



0.308 


0.307 


1.01 • 


10-1 


3.49 


3.47 


1.15 


0.640 


0.647 


9.96- 


10-2 


7.42 


7.44 


1.16 


2.98 


2.83 


9.48- 


10-2 


37.0 


35.7 


1.18 


1.54 


1.39 


9.45- 


10-2 


17.0 


15.7 


1.04 


20000 K 


2.32 • 10-2 


2.32 ■ 10-2 


1.02 ■ 


10-1 


0.273 


0.272 


1.20 


2.25 • 10^2 


2.09 ■ 10^2 


9.83- 


10-2 


0.224 


0.183 


9.79 • 10-2 


0.300 


0.284 


1.01 • 


10-1 


2.81 


2.31 


9.46 • 10-2 


0.358 


0.346 


9.96. 


10-2 


3.41 


2.83 


9.49 • 10-2 
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Fig. 10. — Kx{a)/K.p vs. wavelength. Shown are curves for sihcate grains with different sizes, 
as indicated. Also shown is the black body spectrum for Tes — 2500 and 20000 K. The 
Planck mean opacity was calculated for 1000 K. The unit of A^a is arbitrary. 
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Fig. 11. — (dn/da) K\{a) vs. wavelength. Shown are the curves for sihcate grains with 
different sizes, as indicated. Also shown is the black body spectrum for T^s = 2500 and 
20000 K. The size distribution dn/da is a power-law, with index q = —3.5. The units of 
{dn/da) n\{a) and XBx are arbitrary. 
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Fig. 12. — Size distribution models for case A, with 3] = 1000 K. Shown are the results for 
two envelope models of cosmic silicate grains with a MRN size distribution and size range 
as indicated. The left and right panels show the results for a given stellar temperature. For 
each model, the optical depth was adjusted so that T^ep — 0.1. 
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Fig. 13. — Size distribution models for case B, with Ti = 1000 K. Shown are the results for 
two envelope models of cosmic silicate grains with a MRN size distribution and size range 
as indicated. The left and right panels show the results for a given stellar temperature. For 
each model, the optical depth was adjusted so that Tj-ep = 1- 




Fig. 14. — Extinction optical depth, ta, vs. wavelength for the 2500K models shown in 
Figure 13. Note the similar shape in the optical. 
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Fig. 15. — Size distribution models for case B, with Tj = 1000 K. Shown are three envelope 
models of cosmic silicate grains with a MRN size distribution and size range as indicated. 
The left and right panels show the results for a given stellar temperature. For each model, 
the optical depth was adjusted so that r^ep = 1- 
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Fig. 16. — Extinction optical depth, t\, vs. wavelength for the 2500K models shown in 
Figure 15. 
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Fig. 17. — SEDs for size distribution models. The figure compares the SED for a model with 
a — 0.005 to 0.25 //m with the SED for a model with a — 0.05 to 0.25 jim.. Both models have 
Trep = 0.1 and Ti = 1000 K. 



